Some properties of the resonant state in quantum mechanics and its computation 
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The resonant state of the open quantum system is studied from the viewpoint of the outgoing 
momentum flux. We show that the number of particles is conserved for a resonant state, if we 
use an expanding volume of integration in order to take account of the outgoing momentum flux; 
the number of particles would decay exponentially in a fixed volume of integration. Moreover, 
we introduce new numerical methods of treating the resonant state with the use of the effective 
potential. We first give a numerical method of finding a resonance pole in the complex energy 
plane. The method seeks an energy eigenvalue iteratively. We found that our method leads to 
a super-convergence, the convergence exponential with respect to the iteration step. The present 
method is completely independent of commonly used complex scaling. We also give a numerical 
trick for computing the time evolution of the resonant state in a limited spatial area. Since the wave 
function of the resonant state is diverging away from the scattering potential, it has been previously 
difficult to follow its time evolution numerically in a finite area. 

PACS numbers: 03.65.Yz, 73.22.Dj, 03.65.Nk, 05.60.Gg 



I. INTRODUCTION 

The resonance has been studied extensively for quite a 
long time ^ 0, i i, S i, 0, S i, [M HH, El, [H Q El, P ■ 

It appears in almost every field of physics from classical 
mechanics to quantum mechanics. In spite of this fact, 
however, many fundamental aspects remain to be inves- 
tigated. Particularly in condensed-matter and statistical 
physics, many textbooks indeed introduce the complex 
eigenvalue of the resonant state only phenomenologically. 

Meanwhile, the resonant phenomenon is of increasing 
importance especially in the quantum mechanics of meso- 
scopic devices. When we use nano-devices, we inevitably 
attach leads to them. Hence the devices are always open 
systems and have resonant states; an electron comes into 
the device through a lead, is trapped in the confining 
potential of the device for a while with a finite lifetime, 
and may come out through another lead. This resonant 
conduction has been intensely studied experimentally; for 
example, the Fano resonance [13, El, El, US] has attracted 
much attention. 

In the present paper, we study the resonant state of 
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the open quantum system in two parts. In the first part, 
from Section |ll] to Section IIV) we discuss the physical 
significance of the imaginary part of the resonant eigen- 
value in a somewhat pedagogical manner, having in view 
condensed-matter and statistical physicists as a target 
reader. In the course of presenting our results in various 
places, we have often received many questions regarding 
the well-established knowledge on the resonant state from 
condensed-matter and statistical physicists, and hence we 
find it necessary to review here fundamental properties 
of the resonant state for self-containedness. 

We nevertheless include two new points in the discus- 
sions. First, we show that 

the imaginary part of the energy expectation 
is proportional to a momentum flux going out 
of the system for arbitrary wave functions. 

The relation has appeared in the literature but for each 
resonant state [2l|, [H, , in which case the imaginary 
part of the energy expectation is reduced to the lifetime, 
or the half-width of the resonance and the momentum 
flux is reduced to the real part of the wave number of the 
eigenstate. 

Second, we show that 

the number of particles in a resonant state is 
conserved when we count the number in an 
expanding volume. 
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The resonant state represents a decaying state and its 
particle number decreases exponentially, if we count the 
number in a fixed volume. Since the decay rate is re- 
lated to the outgoing momentum flux as stated above, 
we can keep track of leaking particles by expanding the 
integration volume constantly. As far as we know, the 
above fact has never been stated in the literature albeit 
its simpleness. 

In the second part of the present paper, from Section IVl 
to Section IVI[ we introduce new numerical methods of 
treating the resonant state with the use of the effective 

potential [H, [H, [M 113, M M M M M MMM- 

A logical consequence of the arguments in Sections |TT] 
and mil is that the eigenfunction of the resonant state is 
diverging away from the scattering potential. For this 
reason, treating resonant states numerically can involve 
difficulties. Naive treatment by chopping off the infinite 
space would fail inevitably; the resonant state never ap- 
pears in a closed system and the diverging eigenfunction 
is never represented in a finite space. A conventional way 
is suppressing the divergence by the complex sca ling or 
the complex rotation [23, M, M, M, SH, \Mm 

Hi, EE EE 113, mi • We here use the effective potential, an 
energy-dependent boundary condition, in order to cut off 
the infinite space. We can thereby find resonance poles in 
the complex energy plane, working in a finite space. The 
present method seeks the energy eigenvalue iteratively. 
We numerically demonstrate for a simple model that we 
can indeed find a complex eigenvalue of the resonance 
pole. We found that our method is a rapidly converg- 
ing method as in the KAM theory [IE HEl in nonlinear 
dynamics or in the Newton method of finding the so- 
lution of the nonlinear equation; i. e. the convergence is 
exponential with respect to the iteration step. We can 
also compute the time evolution of the resonant state 
in a finite space. We numerically demonstrate for the 
Friedrichs model that we calculate the dynamics of the 
central part of the diverging resonant eigenfunction. Nu- 
merical calculations for a ladder lattice finding a quasi- 
stable resonant states with very long lifetime will be re- 
ported elsewhere [sTj . 



II. THE IMAGINARY PART OF THE ENERGY 
EXPECTATION 

In the present section, we show that the imaginary 
part of the energy expectation is proportional to the flux 
out of a volume as well as the lifetime of the resonant 
state. Thus the energy is generally complex in an open 
quantum system. The relation has been obtained only 
for each resonant state previously. 

Hereafter, we discuss the one-dimensional case just for 
explanatory simplicity. See Appendix lA II for the corre- 
sponding expressions in the three-dimensional case. Con- 
sider the one-dimensional one-body Hamiltonian 



where 



2m 



2m dx^ 



and V = V{x). 



(2) 



We assume that the potential operator V has a finite 
support 



fipot = {a; I < a; < /}; 



(3) 



that is, we assume that the potential function is reason- 
ably localized around the origin. We also assume that 
the potential operator is Herniitian; in other words, the 
potential function is real: 



V{x)* = V{x). 



(4) 



We will see below that the kinetic energy is not necessar- 
ily a Hermitian operator. 



A. The imaginary part and tlie momentum flux 

Let us check whether the Hamiltonian ([T]) is a Hermi- 
tian operator or not. We can do this by seeing whether 
the expectation value of the Hamiltonian operator is 
real for an arbitrary wave function ^'{x). We define the 
Hamiltonian expectation as follows: 



^p{x)*'Hip{x)dx. 



(5) 



Note here that the volume of the integration is limited 
to a segment 



n = {x \ -L < X < L} 



(6) 



that includes the support of the potential operator V; 
that is, we assume fl D ilpot, ot L > I. The reason for 
the trick of limiting the integration volume will be self- 
evident below. The usual definition of the expectation 
is recovered in the infinite-segment limit after normaliza- 
tion: 



lim 



(7) 



We now compute the imaginary part of the expecta- 
tion. 



2ilm(Mn 



(8) 



n = ic + v, 



(1) 



and see in what case it remains. Owing to the assumption 
that the potential is a real and localized function, the 
expectation of the potential term is real. Only the kinetic 
term therefore contributes to Eq. ([8]) . The expectation of 
the kinetic term is transformed by the partial integration 
as 

{mij)n = f ^{xmx)dx 

= TT- I 1p'{xyip'{x)dx 

^m J_]^ 

-^w.^y^'(^t^-L- (9) 
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The first term of the right-hand side of Eq. ([9]) is real. 
The second term can be complex, however, which yields 
the imaginary part of the Hamiltonian expectation. 

Subtracting from Eq. ^ its complex conjugate, we 
have 

2ilm{^\n\i:)n 
= -^Re(V'(x)*p^(a:)|^^^-V(a;)*p^(x)|^_J(,10) 



The Hamiltonian of an open system is Hermitian only 
when the flux into the system is balanced with the flux 
out of the system. Conversely, when we assume that 
the Hamiltonian of an open system is Hermitian from 
the beginning, the flux must be automatically conserved. 
We see indeed in the textbook example of the tunneling 
phenomenon that the sum of the reflection flux and the 
transmission flux is equal to the incident flux. In fact, 
this physically reasonable phenomenon of the flux conser- 
vation is, from an algebraic point of view, a consequence 
of assuming that the energy variable is real. 



where p is the momentum operator. In light of generaliz- 
ing the expression to higher-dimensional cases, we denote 
the last line as 



(V'IPnlV')ao = V(a;)*j5-0(x)| 



(11) 

where Pn is the normal component of the momentum on 
the surface cJfi. This represents the momentum flux going 
out of the segment f2. Thus we arrive at 



2m 



Re{i>\pn\ip)dn- 



(12) 



This is a very insightful formula in analyzing resonant 
states. 

We have shown this general result for arbitrary wave 
functions. This relation, however, has been already de- 
rived for individual resonant states; see e.g. Eq. (2.13) of 
Ref. [2l|, Eq. (24) of Ref. [H and Eq. (3.19) of Ref. 
For each resonant state tpres,n, the energy expectation is 
reduced to the complex eigenvalue 



inHoo (V'IV')o 2 



(13) 



as will be reviewed below in Section [IV A I The momen- 
tum flux, on the other hand, is reduced to 



B. The imaginary part and the lifetime 

Now we mention the familiar relation between the 
imaginary part of an eigenvalue and the decay rate of 
the resonant state (see e.g. Ref. [13]). We now consider 
the time-dependent Schrodinger equation 



d 



(16) 



We compute the time dependence of the particle number 
in the segment J7, given by 



Its time derivative is transformed as 



(17) 



L 

-L 



dt dt 

^{x,tYn'^{x,t) - -^{x.t) {ii^{x,t) 



dx 
(18) 



Re{ip\p^\ip)dn- = hk„ 



(14) 



again as will be reviewed in Section llV Al Taking account 
of the normalization factor 1/^2k„ of the resonant wave 
function, we see that the relation (fT2ll is reduced to 



r 

4k„ 



2m 



(15) 



which has been obtained previously [U [S^, and will 
be derived in a different context in Section UlI Bl 

The formula ([12]) claims that the imaginary part of the 
energy expectation indicates a momentum flux going out 
of the system. This leads to the following statement: 

The Hamiltonian (the kinetic term, in partic- 
ular) is a Hermitian operator if the system is 
closed, but is generally a non-Hermitian op- 
erator if it is open. 



We notice that the last line of Eq. is related to 
the momentum flux in Eq. (|12[) as 



dt 



Nn{t) - -Im(*|7i:|*>o 



1 



Re(*|pn|*)ao. (19) 



We thus conclude that the imaginary part of the Hamil- 
tonian expectation is related to two quantities, the mo- 
mentum flux going out of the system and the decay rate 
of the particle number in the system. This is a very nat- 
ural conclusion; the decrease of the particle number in 
the system is due to the leak from the system (Fig. [1]). 

Equation (fT9|l also indicates that dealing with the com- 
plex energy eigenvalue is closely related to dealing with 
the complex momentum of the eigenstate. In fact, we 
have found that in some cases it is more convenient to 
solve the problem in the complex momentum plane than 
in the complex energy plane, which we will demonstrate 
below. 
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FIG. 1: The decrease of the particle number in the system is 
equal to the leak from the system. 

FIG. 2: The S matrix of a scattering potential is defined from 
two scattering solutions as in (a) and (b). The singularities 
III. RESONANT STATE AS A SOLUTION OF of the S matrix occur when the incident waves are missing. 
THE STATIONARY SCHRODINGER EQUATION 



The facts established in the previous section indicate 
that the open quantum system can possess an eigen- 
statc with a complex eigenvalue. We generally refer to 
such a state as a resonant state. In the present section, 
we review that the resonant state is an eigenfunction of 
the Hamiltonian of the open quantum system under the 
boundary condition that we have only outgoing waves. 
The statements in this section are known except for a 
comment on a resonance peak of the transmission prob- 
ability presented in Sections IIII Bl and IIII CI Readers fa- 
miliar with the resonant state may skip to Section [ill Bl 

The perhaps new point in the present section is the 
observation that a resonance peak of the transmission 
probability does not necessarily correspond to a resonant 
eigenstate. We show in a tutorial example, specifically in 
Fig. 4(d) and (e), that the lowest peak of the transmission 
probability, or the conductance, cannot be called a reso- 
nance peak in the sense that the corresponding resonant 
eigenstate is missing. 



A. Defining the resonant state as an eigenfunction 

The quantum resonant state is, in most textbooks, de- 
fined as a singularity of the S matrix. We show that 
this definition is actually equivalent to solving the sta- 
tionary Schrodinger equation under a certain boundary 
condition. Specifically, the boundary condition is to have 
outgoing waves only, no incoming waves. The momen- 
tum is obviously not conserved and hence the eigenvalue 
is complex as is indicated by Eq. (fT2|) . 

We again discuss the one-dimensional case. As in the 
previous section, we have a potential V{x) in a finite 
region 17pot defined in Eq. dS]). Outside the region, we 
have only the kinetic term in the Hamiltonian, 

Vix) = for xiflpot, (20) 
and hence any eigenfunction is composed of exjp{±iKx) 



with 



K 



\/2mE 



(21) 



The argument in the previous section shows that a 
momentum flux going out of the central area (fi defined 
in Eq. ([6])) may yield a state with a complex eigenvalue. 
This motivates us to seek an eigenstate with boundary 
conditions with outgoing waves only: 



V'ros(a;) 



for X < — L, 
for X > L. 



(22) 



This set of boundary conditions is often called the Siegert 
condition 0, and has been used occasionally in the 
literature. An eigenfunction satisfying Eq. ([22]) . if it ex- 
ists, obviously has a finite momentum flux going out of 
the central area, and hence its energy eigenvalue E must 
be complex because of the formula (|12p . 

In the present paper, we define the resonant state in 
the following way: 

A resonant state is a solution of the 
Schrodinger equation (an eigenfunction of the 
Hamiltonian) under the boundary conditions 
that only the outgoing waves exist outside the 
segment 51. 

We show below that this definition is equivalent to seek- 
ing poles of the S matrix of the potential. 

In order to define the S matrix, we consider an eigen- 
function satisfying the following boundary conditions 
(Fig. ID) instead of Eq. 



■0(2;) 



J^piKx 
(j^iKx 



Be 



-iKx 



for X < —L, 
for X > L. 



(23) 



In order for Eq. ((23|) to be a solution of the stationary 
Schrodinger equation 'Htl;{x) — Eip{x), the coefficients 
A, B and C must depend on E. The momentum flux out 
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of the central area H can vanish for this set of bound- 
ary conditions whenever the outgoing flux is balanced by 
the incoming flux as in \A{E)\^ = \B{E)\^ + \C{E)\^. 
Then, an eigenstate satisfying Eq. has a real energy 
eigenvalue. 

The S matrix in one dimension is defined as follows. 
The reflection amplitude from the left back to the left is 
given by 



B{E) 
ME)" 



(24) 



while transmission amplitude from the left through to 
the right is given by 



tRLiE) 



C{E) 

ME)- 



(25) 



Likewise, we can consider an eigenfunction satisfying 



I „ — iKx 



C'e 
A 



for X < —L, 
+ E'e'^"" for x> L, 



(26) 



and have the reflection amplitude from the right back to 
the right and the transmission amplitude from the right 
through to the left as 

rME)^^ and t^E)^^. (27) 
The S matrix in this case is then defined in the form 



(28) 



which relates the incoming waves from the left and the 
right, to the outgoing waves to the left and the right. 

Many textbooks define resonant states phenomenolog- 
ically as singularities of the S matrix analytically con- 
tinued onto the complex E plane; the singularities are 
not associated to the true complex eigenstates of the 
Schrodinger equation. The elements of the S matrix in 
fact diverges whenever 



A{E) = or A'{E) = 0. 



(29) 



We notice that the set of boundary conditions ([^5]) is re- 
duced to Eq. when A{E) = 0. Instead of solving the 
Schrodinger equation under the boundary conditions (|23p 
and then seeking the zeros of A{E), Eq. ([22]) means that 
we solve the Schrodinger equation with A(E) — from 
the very beginning. Thus we see that the phenomenolog- 
ical definition of the resonant state given in many text- 
books is equivalent to our definition of the resonant state 
as an eigenfunction. 

Some might think that the outgoing waves in Eq. 
spring out of nowhere; it is of course not true. The par- 
ticle as the source of the outgoing waves is there in the 
trapping potential from the very beginning. The par- 
ticle probability leaks from the potential at a rate pro- 
portional to the particle probability itself; notice that the 



right-hand side of Eq. (fTQ]) has the factor exp{—tlmE/h). 
Hence the particle probability in the trapping potential 
decreases, so does the leaking probability. The parti- 
cle probability keeps decreasing exponentially forever, as 
the leak keeps decreasing forever; see Sec. IIVI below for 
more details. Since we here define the resonant state as 
a stationary state, we can only describe the resonance as 
above. In this picture, we do not employ the dynamic 
description of the resonance that a particle comes into 
the scattering potential from outside, is trapped by the 
potential for a while before it escapes. 



B. Tutorial example 

We here confirm the above arguments, solving a simple 
example in one dimension. Let us consider the potential 
with two delta functions: 



V{x) = Vo{S{x + l) + 5{x-l)). 



(30) 



In order to obtain the resonant state for the potential, 
we assume the form 



fj^iKx 



for X < —I, 

for -I <x <l, (31) 

ioi X > I. 



Note that only outgoing waves exist outside the potential. 
The connection conditions at x = ±1 are given by 

iP,,si±l + 0)-Acs{±l-0) - 0, (32) 
V'L(±^ + 0)-V^;es(±^-0) - a-Vrcs(±0, (33) 

where the length scale a is defined by 



(34) 



2mVo 



The solution PT|) has four undetermined constants (the 
energy eigenvalue E and the three ratios of the coeffi- 
cients B, C, F, G), while we have four connection condi- 
tions Hence we can fix the energy eigenvalues 
to discrete values. 

The solution is simplified further. Since the poten- 
tial ([50]) is an even function, the solutions can be classi- 
fied according to their parities. The even and odd solu- 
tions are given by 



B^±C and G = ±F, 



(35) 



respectively. The solutions of a definite parity have two 
undetermined constants (the energy eigenvalue E and the 
ratio between C and F) , while we have now two indepen- 
dent connection conditions at x — +1. Eliminating the 
coefficients C and F, we arrive at 



1 — 2iKa = =Fe 



2iKl 



(36) 



where the upper sign gives the even solutions and the 
lower sign gives the odd ones. This equation determines 
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discrete wave numbers Kn and hence does discrete energy 
eigenvalues = h^K^/2m in the complex energy plane. 
For numerical calculations, we set 

where A:„ > 0. Equation (j36p is then rewritten as 



I 



Vn = TT 



6> 



2a tan 2^„ ' 



2a 

1 - -nn 



(38) 
(39) 



where ^„ = knl and ?7„ = We plot in Fig. [3] the 

functions ((38)) and ([39)1 in the region ^ = fc^ > in the 
cases a/l = 0.1, 1, 4 and 10. The crossing points are the 
solutions; the circles denote the solutions with even par- 
ity and the squares denote the solutions with odd parity. 
There are also solutions in the region ^ — kl < Q with 
the negative real part — ^„ and the same imaginary part 
r}n\ Eqs. (f38|) and (l39|) are symmetric with respect to the 
T] axis if we take into account the sign of the right-hand 
side of Eq. ([39|. 

The lowest even solution disappears for a/l > 
3.59112- ••. The curve of Eq. ^ splits into two 
branches for a/l > 3.59112 one circling near the 
origin, the other continuing to infinity. In the cases 
Fig.[3](d) and (e), there is a small branch near the origin, 
which does not intersect with Eq. ((38|) at all. The rem- 
nant of the resonance peak corresponding to this missing 
pole nevertheless exists in the transmission probability; 
see the last paragraph of Section IIII CI below. We have 
never seen this phenomenon reported previously. 

Notice here that all solutions (the crossing points in 
Fig. [S]) satisfy 



= -Imi^„ > 0. 



(40) 



We list in Table U numerical estimates (by Mathematica) 
of the solutions Kn and the resulting energy eigenvalues 



En 



2m 



.r„ 



(41) 



The solutions with kn = YieKn > have outgoing 
waves only, whereas the solutions with kn < (not shown 
in Fig. [3]) have incoming waves only. We note here that 
Eq. gl]) gives 



2m 
m 



(kl - kI) and 



(42) 
(43) 



The latter equation is equivalent to Eq. (fT5|) . Since we 
always have 



IniKn < 0, 



(44) 



we have hnEn < 0. or 



r„ > for the solutions with fc„ > 0. 



(45) 



This indicates that the states with outgoing waves only 
arc indeed decaying. For the solutions with kn < 0, on 
the other hand, we have ImEn > 0, or 



r„ < for the solutions with fc„ < 0. 



(46) 



This indicates that the states are growing because of the 
incoming waves; such states are often called anti-resonant 
states. We argue these facts more in Sec. IIVI below. 

A decaying resonant state and a growing anti-resonant 
state always appear in a pair. Both of the states break 
the time-reversal symmetry. The latter states are the 
time-reversed states of the former. 



C. Relation to the 5* matrix 

Let us now demonstrate for the double delta poten- 
tial (|5D|) that the pole of the S matrix gives the same 
answer as those in Table HI Computing the 5* matrix, or 
computing the transmission and reflection amplitudes, 
differs from the computation in Sec. IIIIBI in two points: 
First, the wave function assumes the form 



iKx 
iKx 



for X < ~l, 
for —l<x<l, 
for X > I. 



(47) 



rather than Eq. (|3ip . Second, the wave number K (and 
hence the energy eigenvalue E) is a given real number. 
(We later carry out analytic continuation onto the com- 
plex K plane.) The undetermined constants are four ra- 
tios of the coefficients A, B, C, F, and G, while we have 
four connection conditions Hence the four ra- 

tios can be obtained as functions of K. We can rephrase 
this point as follows: we have five undetermined con- 
stants, the wave number K and the four ratios, whereas 
we have only four conditions, and hence the solutions 
exist continuously on the real K axis. 

The connection conditions ([32)) -((33 )) relates the five 
coefficients in the forms 



[iK 



-iKl 



iKe 



iKe''^^ 



-iKl 



-iKe 



-iKl 



^^Kl 

-iKe'^^ 



(48) 



e 

{iK -a 



iKl 



-l^^^Kl 



C. (49) 
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FIG. 3: Plots of Eq. ((38]) (the solid lines) and Eq. ^ (the broken lines) in the cases (a) a/l = 0.1, (b) a/l = 1, (c) a/l = 3.5, 
(d) a/l = 4 and (e) a/l — 10. The vertical thin gray lines indicate kl = nn, which are the bound states at Vb — * cxa or a = 0. 
The circles indicate even solutions and the squares indicate odd solutions. Note that the first even solution disappears in (d) 
and (e). 
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TABLE I: Numerical estimates of the complex wave numbers and the complex eigenvalues of the first few resonant states of 
the double delta potential (|3Up . The solutions with KeK^ > and Imi5„ < are shown in Fig. |31 Note that the first even 
solution disappears in the cases a/l = 4 and a/l — 10. 



a/l state (parity) 


jy- I. 




J^n — En — *i n/^ ['1 / TTll \ 


0.1 


1st (even) 


±1.4309486581029545770 


- iO. 0180132370706695616 


1.023644792708441123 =F iO. 025776017414365005 




2nd (odd) 


±2.8775774584575874315 


- i0.0665106724899688980 


4.138014179934080207 =F i0.191389611903989679 




3rd (even) 


±4.3478216485135269076 


- i0.1331827632067771627 


9.442907719433745119 =F i0.579054901079276581 




4th (odd) 


±5.8413795860760520688 


- iO. 2064800963021565454 


17.03954071922854077 =F il. 206128619470434570 


1 


1st (even) 


±0.8940940206918146011 


- iO. 3025104586463533055 


0.353945770123213981 =F iO. 270472792272442937 




2nd (odd) 


±2.2985790066512866386 


- i0.7660460609931899527 


2.348319441127416761 =F il. 760817393926857498 




3rd (even) 


±3.8592068943854588960 


- il. 0264132410357959781 


6.919976856149325741 =f i3. 961141056293867654 




4th (odd) 


±5.4340030287668464008 


- il. 1969911205792216173 


14.04780058695087972 =f i6. 504453374634511715 


3.5 


1st (even) 


±0.1281226970608689462 


~ iO. 6318653191999576592 


-0.191419178052756592 ^ iO. 080956288875125433 




2nd (odd) 


±2.081 1 197436940274902 


- il.4192306059503943142 


1.158421937363385626 =F i2.953588834898203941 




3rd (even) 


±3.7327941962053519687 


- il. 6702340594138141162 


5.572035348999205673 =F i6.234640003284390347 




4th (odd) 


±5.3444838317948077028 


- il.8348671938103419869 


12.59838490469733564 =f i9.806418050810082608 


4 


1st (odd) 


±2.0634804406374274361 


- il.4929932341806764916 


1.014461365791977689 =F i3.080762336735840215 




2nd (even) 


±3.7223094302121584017 


- il. 7400125856792561725 


5.413971847962076268 =F i6.476865256361736495 




3rd (odd) 


±5.3369637981977019936 


- il. 9033720492900427597 


12.43017871262713230 =F ilO. 15822772156233025 
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1st (odd) 


±1.9643116049421679709 


- i2. 0079089942082502035 


-0.086589223855955588 T i3.944158938991022210 




2nd (even) 


±3.6591545130696776716 


- 12.2219369473965324296 


4.226203976156184490 =F iS. 130410608822284633 




3rd (odd) 


±5.2907672454808363886 


- 12.3749885684814657674 


11.17582367271761742 =f il2. 56551172651315939 



By eliminating the vector 



we have 



-1 



{2iKa-lf - e 



4iKl 



{2iKa~ l)e 



-2iKl 



{2iKa+l)e 



2iKl 



The S matrix is then given by 



S{K) 



r{K) t{K) 
t{K) r{K) 



with 



r{K) 
t{K) 



B _ iiKacos2Kl + 2ism2Kl 
A ^ {2iKa - 1)2 - e4*^' ' 
C _ -4K^a^ 
A ^ {2iKa - 1)2 - e^*-^' ' 



C. (50) 

(51) 

(52) 
(53) 



Now we carry out the analytic continuation onto the 
complex K plane. The poles of the elements of the S 
matrix are determined by the zeros of the denominator, 



{2iKa- 1) 



4iKl 



0, 



(54) 



which is equivalent to Eq. (|36p in Sec. HUB] (It is notable, 
however, that the parity of the eigenstate is not obtained 
from the S matrix.) 

Note in Eq. (|50p that the coefficient A vanishes at the 
complex zeros of Eq. and the wave function ([Tf|) is 
reduced to the resonant wave function ((3T|) . (The flux is 
conserved only when K is real.) Thus we show that the 
solutions in Table [J are identical to the resonant states 
given by the conventional definition. 



The transmission probability T = |t(if)|2 has "reso- 
nance peaks" on the real k axis (and on the real e axis) 
as the skirts of the poles of the function t{K). In Fig. 21 
we superimposed the locations of the resonant states on 
the k dependence of the transmission probability T. A 
resonant state with a small imaginary part yields a reso- 
nance peak near its real part, but a resonant state with 
a large imaginary part and the corresponding resonance 
peak are relatively apart from each other. 

Note in Fig. [H^d) and (e) that the first "resonance 
peak" does not have a corresponding resonant state as 
mentioned in Section IIII Bl perhaps we should not call 
this peak a resonance peak. 



IV. DIVERGING EIGENFUNCTION AND 
PARTICLE-NUMBER CONSERVATION 

Following the above tutorial example, we now discuss 
the resonant wave function in general cases. We first 
show that the resonant eigenfunction is diverging away 
from the scattering potential. In spite of, or rather, 
thanks to the divergence, the resonant wave function in- 
deed conserves the particle number. 



A. Divergence of the resonant eigenfunction 

The general form of the one-dimensional resonant wave 
function (1^^ now reads 



Be ^'^"^ '^"^ for X < —L, 



(55) 



9 



3 
2 
] 

K/ 
1 

-2 

-3 



— 1 










J 




IJ 





0.98 r 




FIG. 4: The fc dependence of the transmission probability T (the right axis of each panel) in the cases (a) a/l — 0.1, (b) a/l — 1, 
(c) a/l — 3.5, (d) a/l — 4, and (e) a/l — 10. The vertical thin gray lines indicate kl = rni. The locations of resonant states 
given in Fig. |3] are superimposed (the left axis of each panel) . (The circles indicate the even solutions and the squares indicate 
the odd solutions.) Note in (d) and (e) that the first "resonance" peak does not have a corresponding resonance pole. 
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We argue in the following that the above eigenfunction 
is diverging away from the potential. 

The solution (|55p is an eigenfunction of the stationary 
Schrodinger equation under the boundary conditions of 
the outgoing waves only: 



^ 2 ) ^rcs,n- 



(56) 



The solution of the time-dependent Schrodinger equa- 
tion (jl6p is then given by 



i{x,t) = Vros,n(a:) exp 



if .r„ 



(57) 



Then the particle number of this particular state in the 
segment ft, 



iV„,o(i) EE (*rcs.n|*rcs,„)n, 

decays as in Eq. or 



(58) 



dt 



_e-r"*A X i^iV„.^(0). 



(59) 



If the eigenfunction has outgoing waves only, or fc„ = 
ReKn > 0, then the state should decay in time, or 
r„ > 0, because the particles leak from the central seg- 
ment ri. If the eigenfunction satisfies fc„ < 0, on the other 
hand, the state should grow in time, or r„ < 0, because 
the particles gather into the central segment fl. As il- 
lustrated in Sec. IIII Bl and particularly in Eqs. ((42l) - (|45| . 
the above physical conditions 



fc„ > ^ r„ > 0, 
fc„ < ^ r„ < 0, 

are necessary and sufficient conditions of 

Kn = ~ Im Kn > 

owing to the exact relation p^ . or P5)) . 

2fP 



(60) 
(61) 



(62) 



(63) 



Upon applying the inequality (j62p to the wave func- 
tion (|55|) . we realize that the resonant eigenfunction must 
be diverging as \x\ — > oo. 

We can also show the necessity of the divergence of 
the resonant eigenfunction (I55p by the following dimen- 
sional analysis. The formula (fT2|) reads for the eigenfunc- 
tion dniD, 



111 kfi 



rcs,n I YYCs.n res,n |'^res,n/ dVL 



2m 



(64) 



As we expand the volume of the integration, fi, the left- 
hand side of Eq. would increase in the order of the 



volume of the region whereas the right-hand side 
would increase in the order of where d is the 

dimensionality of the space. This apparent inconsistency 
comes from the implicit assumption that the wave func- 
tion is almost constant over the volume of the integration. 
A way out of the inconsistency is that the wave function 
in Eq. (j64p grows or decays exponentially as we expand 
the integration volume fJ. 

Indeed, the solution ((55| gives the left-hand side of the 
formula as 



r r 

res,n | v-^res,n 
^ ^ -l-L 



(65) 

where we have neglected the details of the wave function 
around the origin. Meanwhile, the solution ([55)1 gives the 
right-hand side of the formula (|64|) as 



2m 



rcs,n |'V^rcs,n/5f2 



(66) 



which are equal to each other thanks to the relation (|63p . 
(In fact, it is customary to normalize the diverging wave 
function by first introducing a Gaussian convergence fac- 
tor exp(— ax^) in the volume integration in Eq. (j65p and 
then, after integration, putting a to 0. This gives the 
normalization factor 1/-\/2k„ for the eigenfunction.) 

Thus we here characterize the resonant state as an 
eigenfunction in the form of an exponentially divergent 
outgoing wave and an eigenvalue with a negative imagi- 
nary part. We show in the next subsection that the diver- 
gence is physically plausible; the resonant eigenfunction 
must diverge in order to conserve the particle number. 
How to deal with the divergence is also the central issue 
of Sees. |V] and |Vll we propose methods of getting rid of 
the divergence computationally in Sections fVl and IVTl 



B. Particle- number conservation 

We here show that we can intuitively understand the 
divergence of the resonant eigenfunction from the point 
of view of the particle-number conservation. Because the 
energy eigenvalue is complex, the particle number is not 
conserved in the conventional sense. We show that the 
particle number in a volume that expands in time is in- 
deed conserved. We then see that the resonant eigenfunc- 
tion must be diverging in order to conserve the particle 
number. 

As we show in Fig.[2l the number of the particles in the 
trapping potential decreases because some of the particles 
escape the potential. The particles that escaped run at 
the velocity hkn/m. The idea here is to count the particle 
number in a volume that expands at the same velocity; 
in the simple one-dimensional case, for example, let the 
boundary L expand as 



L{t) 



-t. 



m 



(67) 



11 



Hamiltonian by 



L{t) 



FIG. 5: The wave amplitude decreases exponentially at every 
point, whereas it increases exponentially as |a;| — > cxd. The 
integral is constant. 



Then the particle number included in the expanding vol- 
ume Q,{t) should be conserved. 

We can estimate the particle number ([58]) as 



^n,0(t)(i) 



L(t) 



L(t) 



|'0rcs,n| dx 



-r„t/h J_ 2K„L{t) 



e " ' X — e 



— exp 



m 



(68) 



where we have neglected the details of the wave function 
around the origin. This is indeed constant because of 
the relation The exponential temporal decrease of 

the wave amplitude is supplemented by its exponential 
spatial increase (Fig. [5|) . 

In short, the divergence of the resonant eigenfunction 
indicates that the particle eventually runs away to infin- 
ity. The divergence is not just plausible but necessary for 
the particle number to be conserved. 

We note that the edge of the expanding volume is well 
defined and never spread in spite of the fact that we have 
a nonlinear dispersion relation E = h?K^ /2m. This is 
because we are considering here properties of a single 
resonant state, and hence are not considering spreading 
of an arbitrary wave packet. 



V. NUMERICAL METHOD OF FINDING 
RESONANT EIGENVALUES 

As we argued in Section IIVI the eigenfunction of the 
resonant state necessarily diverges away from the scatter- 
ing potential. It can be difficult to treat resonant states 
numerically because of this divergence. The resonant 
state would disappear if we naively cut off the infinite 
space for numerical calculations. A conventional method 
is to suppress the divergence by the com plex scalingjOr 
the conmlex rotation [23, [H, IM IH, [H, EO, El, El Si, 
H m, EE 113, The complex rotation modifies the 



P 

X 



p{e): 

x{9) 



; pe 
: xe 



-iB 



(69) 
(70) 



The modification of the momentum operator, Eq. (|69p . 
makes the resonant eigenf unctions in the L"^ functional 
space. The modification of the coordinate operator, 
Eq. (|70|) . which must be introduced to conserve the com- 
mutation relation, can be problematic, however. The po- 
tential V{x) must be an analytic function in order for 
the modification ([70)1 to be apphed [36, 37, 38]; we can- 
not apply the complex scaling to a box potential nor to 
any lattice models. When the potential V{x) is a Gaus- 
sian scattering potential, the complex scaling is applica- 
ble only for < 7r/4; otherwise, the potential becomes 
divergent. 

In the present section, we propose a new numerical 
method of finding resonant eigenvalues iteratively in the 
complex energy plane. Our method is completely inde- 
pendent of the complex scaling and free from the above 
restrictions; in fact, we demonstrate the method for a lat- 
tice model. The method utilizes the effective potential, 
which is called the self-energy of the leads in the context 
of the condensed-matter physics. The effective potential 
is nothing but an energy-dependent boundary condition 
and is useful in cutting off the infinite space of the open 
system. 

The resonant eigenvalue is exactly obtained once the 
energy-dependent effective potential is fixed exactly. To 
fix the effective potential exactly, however, we need to 
know the resonant eigenvalue exactly. In other words, 
the resonant eigenvalue and the effective potential must 
be self-consistent. We make a postulated value of the 
eigenvalue converge to the exact value by seeking the 
self-consistency between the eigenvalue and the effective 
potential. 



A. Effective potential 

Before presenting our numerical method, let us in- 
troduce the effective potential by following Ref. [sst . 
The effective potential was first introduced in nuclear 
physics half a century ago [H, [H, [23] . It is used in 
physics of mesoscopic systems from time to time in recent 
years jH, [H, US [Sll, Hi [H, H [p . 

We demonstrated in Section IIIII that we can obtain 
the resonant eigenvalue using the boimdary condition of 
purely outgoing waves. When the scattering potential is 
complicated, we may need to solve the Schrodinger equa- 
tion numerically. Numerical solution of the Schrodinger 
equation often begins with discretization of the space, if 
the space is continuous. Regardless of whether the sys- 
tem is the tight-binding model or the continuous model 
after discretization, the Schrodinger equation takes the 



12 



form 

H?/'rcs,n(a;) 

= ^ {ipics.nix - Ax) + 1prcs,n{x + Ax)) 

rcs.n 

(x) 

= S„l/'rcs,n(a;), (71) 

where Aa; denotes either the lattice constant or the grid 
size of the discretization mesh. We dropped in Eq. ((7T|) 
the constant potential term ihV'ics," (2^) that arises from 
the discretization. 

Suppose that we try to solve the Schrodinger equation 
in the finite region —L<x<L, which contains the 
support of the potential; in other words, the potential is 
zero at the boundaries and outside the region: 

V{x) EE for |a;| > L. (72) 

Then, the resonant eigenfunction at the boundaries must 
contain the outgoing wave only: 

V'res,„(x) cx e^^"l"l for \x\ > L. (73) 

By using the form (|73p in the Schrodinger equation in 
the region without the potential, we find that the wave 
number Kn is related to the eigenvalue En as 

En^~tY,cos{KnAx), (74) 

where En and Kn are generally complex numbers; see 
Eqs. (I80l)-(l83|) below for the derivation. 
We here cast Eq. ([75)) into the form 

?/;rcs,n (± {L + Ax)) = e*^"^^Vrcs,„ (±i) (75) 

and substitute it for the term in Eq. ([7T|l outside the 
region in question. We thus have 

(±L) 

= - yVA-os,™ {±L - Ax) - y Vros,n (±i + Ax) 

= -y-0ros,ri (± (i - Ax)) 
res,n 

{±L) 

= Enll^rcsA^), (76) 

where the effective potential is given by 

Kff,„(±i:)^-|e'^"^^ (77) 

Using Eq. ([74]) . this is rewritten in terms of the energy 
eigenvalue as 

V,s,n {±L) ^l^En- i\ftJ'^E^^ , (78) 

where we assumed Re Kn > in order to fix the sign of 
the imaginary part. 



Equation (|78|) is the form often found in the litera- 
ture {e.g. Refs. [H, (s^l), where the effective potential is 
obtained as the self-energy that comes into the energy de- 
nominator of the Green's function. The easier derivation 
above is noted in Ref. [s^. 

The effective potential ([78]) is an energy-dependent 
complex potential; the potential is different for differ- 
ent resonant states. As a result, the Hamiltonian matrix 
(such as Eq. ([99]) below) with the effective potential is a 
non-Hermitian operator and may have a complex eigen- 
value. The point is that Eq. ([75]) contains only the eigen- 
functions inside the region —L<x<L. Hence we can 
solve the Schrodinger equation by restricting ourselves to 
the region in question. 



B. Self-consistent solution 

We now propose a numerical method of finding a reso- 
nant eigenvalue by utilizing the self-consistency between 
the resonant eigenvalue and the effective potential. The 
method starts with a postulate of a complex resonant 
eigenvalue E'^^\ This gives the corresponding complex 
effective potential 

V}^^ i±L) = \ [e^') . (79) 

We can then diagonalize the Hamiltonian in the region 
—L < X < L with the effective potential V^g\ From 
the resulting eigenvalues we choose the one closest to the 
postulate i?'"' and refer to it as an updated postulate 
E'-^I This again gives the corresponding effective poten- 
tial V^g\ We repeat the procedure of 

1. obtaining an updated postulate E^'^'> from the di- 
agonalization of the Hamiltonian with the effective 
potential V^^ 

2. and setting the effective potential V^g from the 
postulate 

For a few simple examples that we have tested, the 
convergence is surprisingly good even with the above 
straightforward iteration. A more sophisticated itera- 
tion may be necessary in more complicated cases. In 
the following, we demonstrate that the straightforward 
iteration for a simple model Hamiltonian yields an expo- 
nentially rapid convergence. 

We solve the Schrodinger equation ([7T|l with the po- 
tential 

V{x)^l^'>^ forx = ±Ax, 
^ ' 10 otherwise. ^ ' 

Let us first solve the problem exactly. As we showed in 
Section [nil we assume that the eigenfunction has outgo- 
ing waves only. We have solutions with even parity and 
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solutions with odd parity. Equations ([22)) and (p5l) show 
that we can assume the form 

( ±B„e-'-^"=" for X < -Ax, 

V'res,n(a;) = < Fn for X = 0, (81) 

[ B„e*^"^ for X > Ax, 

where i?„ and F„ are constant. The Schrodinger equation 
for |x| > Ax, 

- y ('0rcs,n(2; " Ax) + ^/-ros.nla: + Ax)) = Enlprcs.n{x) , 

(82) 

gives the dispersion relation 

En = -th cos(ii:„Ax). (83) 
The Schrodinger equation at x = Ax and at x = reads 



2 



(0) + V'res,n(2Ax)) + V {l)lprcs,n{^x) 

(Ax), (84) 

(V' rcs,n 

(-Ax) + ?/'i.os.n(Ax)) 

(0). (85) 

The two equations are reduced to 



for even parity, (87) 

= -F„ (e*^"^- + e-*^"'^-) 

for odd parity, (88) 

where Vq = Vo/th, or 

(2Vb^„ + 1)B„ - F„ - 0, (89) 
2znBn — {zn + Zn^^)Fn — for evcu parity, (90) 
F„ = for odd parity, (91) 

where z„ = e'^"^^. For the odd parity, we have the 
solution 

1 



21^0 



(92) 



1 



if„ Ax = TT + i log 2Vo and £'„ ^ Vq + (93) 

which is a bound state with a real energy. For the even 
parity, we eliminate the degree of freedom of Bn/ Fn, ar- 
riving at the third-order equation 



2VoZn^ - Zn^ + 2%Zn + 1=0. 

The three solutions for Vg = 1, for example, are 

Enlh^ = 1.517526485679543..., 
Enlh^ = -0.383763242839771... 

-i0.132164836187054... 
Enlt\^ = -0.383763242839771... 

-|-i0.132164836187054... 



(94) 

(95) 
(96) 
(97) 




FIG. 6: A three-dimensional plot of log \D[E) \ in the complex 
energy plane. 



The first solution is another bound state. The second 
solution is a resonant state and the third solution an 
anti-resonant state. 

These are the exact solutions of the real and complex 
energy eigenvalues. We now demonstrate our numerical 
iterative method with the rapid convergence to obtain 
the resonant state with the complex eigenvalue. 

We set L — 2 for the numerical procedure proposed 
here. In order to set the first postulate, it may be useful 
to scan the quantity 



D{E) = |det {n,s{E)-EI)\ 



(98) 



in the complex E plane, where the effective Hamiltonian 
Ti.cs{E) is defined on the region —L<x<L with the 
effective potential; for L = 2, for example, we have 



'Hcs{E) — 



( V,s{E) 


2 
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th 
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th 

2 


V,s{E) ) 



(99) 



where 



Kff(i?)=-(i?-zW-i?2 



(100) 



The zeros of D{E) indicate the eigenvalues, although it is 
not practical to use the zeros as numerical estimators of 
the eigenvalues. Figure [5] shows a three-dimensional plot 
of log|D(£')| in the complex energy plane. The dimple 
in the lower left corner is the resonance pole ([96|l . 

On the basis of rough estimation in Fig. [6l we set the 
first postulate as E'-°^ = -0.3 - iO.l. The postulate 
converged to the eigenvalue ([96]) exponentially with re- 
spect to the number of iterations, q; see Fig. [71 The 
numerical calculation involves diagonalization of 5 x 5 
non-Hermitian matrices. 
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FIG. 7: The convergence of the numerical procedure proposed 
here for the model ((80l) with Vb/th = 1. The postulate 
approaches the exact value exponentially with respect to the 
step number q. 

VI. NUMERICAL COMPUTATION OF 
DYNAMICS OF RESONANT EIGENFUNCTIONS 

We finally introduce a trick for numerical computation 
of the dynamics of diverging resonant eigenfunctions with 
the use of the effective potential. The trick introduced 
below enables us to calculate the time evolution of the 
divergent eigenfunction precisely in a finite space, when 
we discretize the space for the numerical calculation or 
when we consider a tight-binding model. 

Now, we are interested in integrating the time- 
dependent Schrodinger equation 

d 

in— *res,«(a;, t) = n^,^^^n{x, t). (101) 

In the right-hand side, we again cutoff the space |a:| > L 
with the use of the effective potential: 

+ VeS,n {±L) -frcs.n {±L, t) . 

(102) 

Since we already know energy eigenvalues by the method 
in Section|Vl we can fix the effective potential Vcs^n {^L) 
for each resonant state. The point is that Eq. p02p con- 
tains only the eigenfunctions inside the region —L < x < 
L. Hence we can integrate Eq. (jlOip by restricting our- 
selves to the region in question. 

Let us demonstrate the above trick in the computation 
of a version of the Friedrichs-Fano (or Newns- Anderson) 
model (Fig. [H]) [H, m. Hi, [13 : 

n - -|^(|x + Ax)(.T| + |a:)(a; + Aa:|) 

X 

-H5(|d)(0K|0)(d|)-|-£;d|d)(d|, (103) 

where d denotes the site of an adatom. The first term 
denotes the tight-binding hopping on a chain, the second 



d 



-3-2-10 1 2 3 
FIG. 8: The tight-binding model on a chain with an adatom. 

term denotes the hopping between the origin of the chain 
a; = and the adatom, and the third term denotes the 
impurity level of the adatom. 

Let us first find the resonance poles exactly. As we 
showed in Section IIIIl we assume that the eigenfunction 
has outgoing waves only. We have solutions with even 
parity and solutions with odd parity. The odd solutions 
vanish at x = 0, and hence do not couple to the adatom. 
We consider only even solutions hereafter. We assume 
the form 

AesA^) = ^ne^"^"'"' (104) 

for the chain and 

'^res,n 

(d) = F„ (105) 

for the adatom. The same procedure as in Section fVB I 
yields the fourth-order equation 

+ 2EdzJ^ + 4g2z„2 _ 2EdZn -1 = 0, (106) 

where z„ — e*^""^^, = E^/th and g = g/t^- Hence 
there are four solutions, as shown in Fig. [9] for g — 0.1. 

Through the dispersion relation En — — cos(if„Aa;), 
or the relations 

th(^^„ + ^^ = -2E.n, (107) 

th (^zn - ^) = ±2^E„^-th^, (108) 

ih' (^zn^ - = t4S„/e7^V, (109) 

the fourth-order equation (|106p is transformed to an 
equation with respect to the energy: 

± {En - Ed) ^fsj^t^ ^ g\ (110) 

This is equivalent to the dispersion equation obtained in 
Ref. [13] for the Friedrichs model; see the pole of Eq. (7) 
of the article. (The correspondence of the notations is 
as follows: En <^ z, Ed ^ Eq, ti^ 4^ B, g ^ gB.) 
We stress that the equation is easier to obtain in the 
momentum space than in the energy space, because we 
need not perform a complicated integration of the self- 
energy part as in the case of the energy space. 

Two of the solutions (the solid and broken lines, a and 
b, in Fig. (5]) are bound states (positive Im Kn and real 
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FIG. 9: The four solutions of Eq. (pL06|l for g = 0.1. The de- 
pendence on the impurity level of, from above to the bottom, 
the real part and the imaginary part of {Kn} and the real 
part and the imaginary part of {En\ for g = 0.1. The solid 
line (a) and the broken line (b) indicate bound states, while 
the dotted line (c) and the chained line (d) indicate resonant 
states. 
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FIG. 10: The four solutions for g = 0.1 and £'d = -1/2 on the 
complex momentum plane and the complex energy plane. In 
each panel, the symbols a, b, c and d indicate the solutions of 
the solid line, the broken line, the dotted line and the chained 
line in Fig. (9] respectively. The solutions a and b are bound 
states, while the solutions c and d are resonant states. 



En) all through the rang e and are referred to as super- 
stable states in Ref. [531 ■ The other two solutions (the 
dotted and chained lines, c and d, in Fig. ^ have nega- 
tive Imifn, which indicates that they are resonant states. 
In some regions, the solutions are pure imaginary and 
hence the energy eigenvalues are real. The states with 
pure irnaginary wave number are often called anti-bound 
states [53l] . but we here include them in resonant states. 

Figure [TU] shows the four solutions for 5 = 0.1 and 
Ed — —1/2 on the complex momentum plane and the 
complex energy plane. The two solutions a and b are 
bound states. The solution a has Reif„ = tt and the 
solution b has YleKn = 0. They are both on the first 
(so-called physical) Riemann sheet of the complex en- 
ergy plane. The other two solutions are resonant states. 
The solution c has JieKn > and ImEn < 0; this is a 
decaying state with an outgoing wave only. The solution 
d has Re Kn < and Im En > 0; this is a growing state 
with an incoming wave only. These two are conjugate 
to each other and always appear in a pair, as was also 
demonstrated in Sec. IIII Bl They are both on the sec- 
ond (so-called unphysical) Riemann sheet of the complex 
energy plane. 

We stress here again that it is generally more conve- 
nient to solve the problem in the complex momentum 
plane than in the complex energy plane. This is because 
of the added complexity in identifying the Riemann sheet 
for a given solution of the complex energy eigenvalue in 
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the case where there are several Riemann sheets. It is 
generaUy not easy to tell whether a state is on the first 
or the second Riemann sheet before we find the location 
of the state in the complex momentum plane. 

The eigenfunctions {ipics.ni^)} of the solutions a, b, c 
and d for g = 0.1 and Ed = —1/2 are shown in Fig. [TTl 
where the normalization is always fixed by '>pTcs,nid) ~ 
1. The eigenfunctions from the top to the bottom of 
Fig.[TT]correspond to the solutions a, b, c and d in Fig.fTUl 
respectively. It is evident that the solutions a and b are 
bound states decaying in space whereas the solutions c 
and d are resonant states diverging in space. 

We then solved the time-dependent Schrodinger equa- 
tion 

d 

ih—'^^cs,n{x, t) = 7Ycff*rcs,n(a;, t) 

-^*rcs,n(2; + ^X,t) 

+ Vc(i.n'^rcs,n{x, t) for X = -L 

(*rcs,n(a; - Ax, t) + *„s,n(a; + Ax, t)) 

+V{x)^!^cs.7i{x, t) for -L < x < L 



2 



(111) 



+V'off,n*rcs,n(a;, t) for X = L 



in the range —L < x < L by the Runge-Kutta method, 
taking the eigenfunctions {V'rcs,n(a^)} as the initial con- 
ditions. The time evolution is partly shown in Fig. [12] 
for g = 0.1 and Ed = —1/2. The bound states a and 
b (the first and the second rows in Fig. [T2)) are time- 
independent and stable, while the resonant state c in the 
third row decays in time and the state d in the fourth 
row grows exponentially in time. 

It should be emphasized that if we would solve the 
original Schrodinger equation simply by truncation of the 
domain of x without introducing the effective potential, 
we could not obtain such accurate numerical results as in 
Figs. [TTjandfT^ This is especially true for the resonant 
eigenstates with the complex eigenvalues, which diverge 
in X. 

We have applied the above technique to a more com- 
plicated model on a ladder lattice. As a result, we have 
found a quasi-stable resonant state with a very long life- 
time for quite a wide range or parameters. The state 
appears to be a bound state around the impurity, diverg- 
ing only far away from the scattering center, and barely 
decays for quite a long time. This finding will be reported 
elsewhere ISll. 
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FIG. 11; The eigenfunctions at t = for g = 0.1 and Ed = 
— 1/2: each panel shows the solution of, from top to bottom, 
the solid line, the broken line, the dotted line and the chained 
line in Fig. [O] or the solutions a, b, c and d in Fig. 1101 The 
solutions a and b are bound states, while the solutions c and 
d are resonant states. In each panel, the absolute value (the 
solid line), the real part (the broken line) and the imaginary 
part (the chained line) are displayed in the region — 20 < a; < 
20. (In the first and second panels, the imaginary part is 
zero.) 



VII. SUMMARY 

In the present paper, we have presented some prop- 
erties on resonant states of open quantum systems in 
two parts. In the first part, we mostly gave a review on 



properties of quantum-mechanical resonant states. We 
nevertheless included two new points: (i) the imaginary 
part of the energy expectation is proportional to a mo- 
mentum flux going out of the system for arbitrary wave 
functions; (ii) the number of particles in a resonant state 
is conserved when we count the number in an expanding 
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FIG. 12: The time evolution of the eigenfunctions at t = 0,20,40,60 for g = 0.1 and Ed = —1/2 is shown in the region 
—20 < X < 20: from top to bottom, the time evolution of the solid line (a), the broken line (b), the dotted line (c) and the 
chained line (d) in Fig.|9] or the solutions a, b, c and d in Fig. IIOI respectively; from left to right t = G, 20, 40, 60. The solutions 
a and b (the first and second rows) are bound states, and hence there is no time evolution in fact. The solutions c and d (the 
third and fourth rows) are resonant states. In each panel, the dot indicates '^rcs,n{d,t) and the curve indicates \'^ies,n{x,t)\. 



volume. We also give an example where a "resonance 
peak" of the conductance is not accompanied by a reso- 
nant state. 

In the second part, we gave numerical methods of find- 
ing a resonant state and computing its dynamics, using 
the effective Hamiltonian. The method of finding a reso- 
nance pole is independent of the complex scaling and is 
applicable to singular potentials and lattice models. This 
method is very efficient, showing the rapid convergence to 
the exact solution. The method of computing the time 
evolution of the divergent resonant eigenfunctions was 
demonstrated for the Friedrichs-Fano model. 

We have intentionally avoided discussing the complete- 
ness relation and the resonant-state expansion [58, 59, 60, 
[sil . [6^ . [63l . [6^ [65j . We are planning to report discussions 
on the topic elsewhere. 
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APPENDIX A: EXPRESSIONS IN THE 
THREE-DIMENSIONAL CASE 

In the main text, we focused on the one-dimensional 
case ([T]). In the present Appendix, we give the cor- 
responding expressions in the three-dimensional case, 
■H = IC + V, where 



/C = — = V^ and V = V(x). 



(Al) 



1. The imaginary part of the resonant eigenvalue 

Instead of Eq. ([5]), we define the Hamiltonian expecta- 
tion as 



'ijj{xyn'iij{x)dv, (A2) 



where we consider an integration volume that includes 
the support of the potential support, fipot- The transfor- 
mation in Eq. ([9]) then corresponds to Gauss's theorem: 

= ^jjj^^^{sr-^^{x)dv 

-|^j^^^(a')*VV'(x)-dS'. (A3) 

In the second term of the last line, the integration is 
carried out over the surface 9il of the volume 17 and dS is 
a vector normal to the surface and of the magnitude of the 
surface element; hence only the normal component of the 
differentiation V contributes to the integral. Subtracting 
from Eq. (|A3ll its complex conjugate, we have 

2zlm(7/;|7i:|i/;)f2 = -— Re // tlj{x)*^j{x) ■ dS, (A4) 
m JJan 



which is followed by the formula (|T2|) . where we define 
(V'|Pn|^)9o= // ^jixYp^ix) ■ dS. (A5) 

The expressions in Sec. IIIBI do not differ much in the 
three-dimensional case. 



2. The S matrix in three dimensions 

We here show in three dimensions the equivalence 
stated in Sec. lIII Al that is, the equivalence of seeking the 
singularities of the 5* matrix and solving the Schrodinger 
equation under the boundary condition of the outgoing 
wave only. 

Let us first review the theory of the S matrix in three 
dimensions. We show that singularities of the S matrix 
can have large contributions to the scattering cross sec- 
tion. Consider scattering on the basis of the Schrodinger 
equation 



.^y^ + Vix)\^ix)^^m- (A6) 



By assuming the solution of the form 

iPix) ~ 6^'=^ + —f(0; k) as If I ^ oo, (A7) 
r 

we have the differential cross section as 

'ii = \mk)\\ (AS) 

Thus we need to obtain the function f{0;k). For the 
purpose, we expand the wave function (jA7[) in the partial 
waves. First, the incident wave is expanded as 

oo 

^^kz ^ ^{21 + j I {kr) Pi {cos 9) (A9) 

-. oc 

— ^(2/ + 1) [e*'^'- - i-lYe-^'^-] PiicosO) 



2ikr 



as r ^ oo, 



(AlO) 



where ji denotes the spherical Bessel function and Pi is 
the Legendre polynomial. We used an asymptotic form 
of the spherical Bessel function in moving from Eq. (|A9p 
to Eq. (jA10[) . Next, the scattered wave is expanded as 



ikr °° 



— fid; k) = -— + l)aiik)Pi{cos9) (All) 



1=0 



with some coefficients {ai{k)}. By summing up 
Eqs. (|A10|1 and (jAlip . we have 



tpix) 



2ikr 



^(2Z + 1) [5z(fc)e''='- - (--l)'e"*''] Piicos9) 

1=0 

as r ^ cx), (A12) 



where 



Siik) = 1 - aiik) 



(A13) 
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is called the S matrix. The cross section 
lated from the S matrix as 



da 



— J2i2l + l)\l-SdkW 



HI) is calcu- 



(A14) 



1=0 



The singularities of the S matrix thus can affect the cross 
section largely. In the conventional theory, the resonance 
is defined as the singularity of the S matrix. 

We next review how we can compute the S matrix. For 
the purpose, we introduce the Jost solutions and the Jost 
functions. We show that the singularities of the S matrix 
is caused by the zeros of a Jost function. Consider a 
solution of the Schrodinger equation with the momentum 
k and the angular momentum I: 



(A15) 



For simplicity, we drop the subscript I for various quan- 
tities hereafter. The equation for x(r; k) is given by 



X{r;k) 



2m^ . , l{l + 



1) 



X(r;fc) = k'^x{r;k), 



(A16) 

where we consider only bounded potentials satisfying 



lim r V{r) ~ 0, and lim rV{r) 

r — ^0 r — ^oo 



0. (A17) 

Near the origin r ~ 0, Eq. (jAl6|) is reduced to 

1(1 



X{r;k) 



1 



X(r;fc), 



(A18) 



and hence there are two solutions of the forms r +^ and 
r^'. The physical solution must be regular at the origin. 
We thus choose the solution which satisfies 



X(r;fc) 



0. 



(A19) 



Far away from the origin, on the other hand, there are 
two potential- free solutions of the forms exp(±i/cr). Let 
f±(r; k) denote the solutions that satisfy 



f±{r;k) ~ e 



ztikr 



(A20) 



These are called the Jost solutions. The physical solution 
x(r; k) is a superposition of these two Jost solutions: 



X(r; k) - a+ik)f+ir; k) + a^{k)f^{r- k) 



(A21) 



with some coefficients a±(k). For later convenience, we 
change the notation as follows: 



X(r;fc) 



1 

^ik 



(/_(fc)/+(r;fc)-/+(fc)/_(r;fc)); (A22) 



that is, 

f^{k) = ±2ifca±(fc), (A23) 
which are called the Jost functions. Note that the Jost 
solutions (jA20[) . near the origin, are in return superposi- 
tions of the solutions of the forms r'+^ and r^': 



/±(r; k) ~ b±{k)r 



- c±{k)r 



(A24) 



with some coefficients b±{k) and c±{k). The superposi- 
tion (|A22p . in the limit r ^ oo, behaves as 



^ I asr— >oo. 

2tk \f+{k) 

(A25) 

Comparing the asymptotic forms (jA25P and (jA12[) . we 
know that the S matrix is calculated from the Jost func- 
tions as 



Siik) = (-1) 



i f-jk) 



(A26) 



Therefore large contributions to the cross section come 
from the zeros of the Jost function f+{k). 

We can show that the Jost functions, or the superposi- 
tion coefficients f±{k) in Eq. (|A22p . is obtained from the 
intercept of the Jost solutions at r = as 

/±(fc) = {21 + l)c±(fc) = {21 + 1) lim r'/±(r; fc). (A27) 

For the present purpose, however, we simply put f+{k) — 
in Eq. (|A22p . which is reduced to 



X{r;k) ~ —f_{k)f+{r:k) cx e' 



(A28) 



This means that we should seek a particular form of the 
solution xl?'; k) under the two boundary conditions (|A19p 
and (jA28p . This is nothing but solving the Schrodinger 
equation under the boundary condition of the outgoing 
wave only. The only difference is that we normally solve 
the wave function x(''; k) outward (starting from the in- 
ner boundary condition (|A19p and ending with the outer 
boundary condition (IA22|) ). while we normally obtain the 
Jost function f±{r;k) inward (starting from the outer 
boundary condition (jA20|) and ending with the inner 
boundary condition (|A24p ). 

Note that in solving the above problem with the two 
conditions (IA19P and (|A28p . we have two undetermined 
variables k and f-{k). Hence we obtain the solutions for 
discrete values of k. 
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